Natural noise and external wake field seeding in a proton-driven plasma accelerator 
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We discuss the level of natural shot noise in a proton bunch-driven plasma accelerator. The 
required seeding for the plasma wake field must be larger than the cumulative shot noise. This is 
the necessary condition for the axial symmetry of the generated wake and the acceleration quality. 
We develop an analytical theory of the noise field and compare it with multi-dimensional simulations. 
It appears that the natural noise wake field generated in plasma by the available at CERN super- 
protons-synchrotron (SPS) bunches is very low, at the level of a few 10 kV/m. This fortunate fact 
eases the requirements on the seed. Our three dimensional simulations show that even a few tens 
MeV electron bunch precursor of a very moderate intensity is sufficient to seed the proton bunch 
self-modulation in plasma. 

PACS numbers: 41.75.Lx, 52.35. Qz, 52.40.Mj 



Particle beam-driven plasma wakefield acceleration 
(PWFA) is capable of producing accelerating gradients 
far in excess of those in conventional accelerators [![, 
but needs the drive beam to be properly shaped or com- 
pressed (see, e.g., ref. for PWFA basics). So far PWFA 
has been experimentally studied with electron or positron 
beams shaped before the plasma [l], lil-fioj. Recently a 
new approach was proposed pr| - [l3j which assumes beam 
shaping by the plasma itself as a result of the transverse 
two-stream beam-plasma instability [l4l - fl7| . At the non- 
linear stage, the instability splits the initially long beam 
into short bunches spaced exactly one plasma wavelength 
apart [ll| [l9j]. Harnessing the instability would make 
it possible to excite strong wakefields by initially long 
beams without a complicated and expensive compressor 
or chopper. The controlled instability is the key physical 
effect to be demonstrated by the discussed proton-driven 
PWFA experiment in CERN [H and auxiliary ex- 
periments [Him. 

To be useful for acceleration of a witness beam, the 
generated wake must be axisymmetric. Yet, the plasma 
supports various modes of the instability, including non- 
axisymmetric (or hosing) ones. The latter are undesir- 
able since they quickly destroy the beam to the state at 
which no strong wakefield is excited (l2| . Fortunately, if 
the axisymmetric mode has grown up to sufficiently high 
amplitude, it prevents development of other modes [ll| . 
Simulations show that the proper mode may dominate 
when an externally introduced seed perturbation is in- 
troduced 12[ . The amplitude of the perturbation must 
be much higher than the noise level from which hosing 
modes grow up. To simulate the instability correctly and 
to determine the required amplitude of the seed pertur- 
bation, we have to know the noise level. 

There could be various sources of uncontrollable seed 
perturbations for instabilities. The one we are interested 



in is the shot noise of individual beam particles. The 
wakefield pattern in this case moves with the beam and 
will be amplified by the beam instability. 

The shot noise field is a sum of wakefields left behind 
by separate beam particles. The contribution of a single 
proton located at the radius r b into the on-axis wake- 
field at the distance z b downstream can be taken from 
Ref. [H: 



E bz 



2ek p K {k p r b ) cos(k p z b ), 



(1) 



\E bJ _\ = 2ek 2 p K 1 {k p r b )&va(k p z b ), (2) 

where k p = u) p /c is the plasma wavenumber determined 
by the plasma frequency u p and the light velocity c, e > 
is the elementary charge, Kq and A"i are the modified 
Bessel functions of the second kind. We use the cylin- 
drical coordinates with the z-axis being the direction of 
beam propagation. By writing the transverse field in the 
form ((2|) we take into account that the wave amplitude 
is low and thus there is no magnetic field left behind the 
particle. 

We will illustrate the obtained formulae by beam and 
plasma parameters discussed in Ref. [19( as SPS-LHC 
variant: number of beam particles N = 1.15 x 10 11 , 
plasma density n p = 7 x 10 14 cm~ 3 , radius ay = fc" 1 = 
0.2 mm, length a z = 12 mm. 

First we calculate the longitudinal noise field E nz on 
the axis. As the beam evolves slowly as compared to the 
plasma timescale, the field depends on longitudinal co- 
ordinate and time in the combination £ = z — ct. We 
consider the field as the real part of some complex func- 
tion: 

EnziQ = Re (£ c *(£)e iM ) ■ (3) 

Each beam particle makes a contribution Qj into the 
complex amplitude E cz which depends on particle loca- 
tion with respect to the observation point. The absolute 
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FIG. 1. (a) Excitation of the noise field as a random walk; 
(b) second moments of probability distributions. 



value of the contribution depends on r^, while the ar- 
gument is determined by z\,. The process of noise field 
excitation thus can be considered as a two-dimensional 
random walk {Rn} on the complex amplitude plane 
(Fig. R N = Q 1 + ... + Q N . 

Properties of the random walk are determined by prob- 
ability distributions of absolute value and direction of 
the steps. Since the beams of interest are much longer 
than the plasma wavelength 2nk~ l , we can assume the 
isotropic distribution of the steps. It is this assumption 
that rules out the collective wakefield excitation and re- 
tains only the shot noise. The distribution of absolute 
values is characterized by the probability density /(s) 
that, in turn, is determined by the radial distribution of 
beam particles: 



at as 



where r(s) is the function reciprocal to 
s(r) = 2ek 2 K {k p r). 



(4) 



(5) 



Thus we have to determine the expectation E az = E|i?/v| 
for large N (the boldface E denotes the expectation). 

Suppose that {Bn} is a sequence of positive constants, 
and the distribution of Rn /Bn converges weakly to the 
two-dimensional normal distribution with the density 
function 



d{ux, v 2 ) = exp I - ~(z/? + v\) 



(0) 



Let £2) be a random vector with the density function 
d(fi, v-i). Then its components are independent and have 
standard normal distribution. Denote p 2 = (f + £|. A 
straightforward calculation shows that p has an expo- 
nential distribution with the density function 



e(" ) = 2 ex P 



v > 0. 



(7) 



We have, therefore, for large N 



E 



Rn 



B N 



Ep = EJC 1+ C. 



ly e(y)dv 



(8) 



Thus, the problem reduces to finding a sequence Bn that 
provides the convergence to the two-dimensional normal 
distribution. 

Let us choose an arbitrary straight line containing the 
origin and denote by <pj the angle between this direction 
and Qj. We suppose that random variables ipj, j = 
1,2,... have uniform distribution on [0, 2tt]. Denote 



\Qj 



Si COS ipi 



■ Sn COS IfN- 



(9) 



We come to the one-dimensional random walk {S'at} 
which is the projection of the initial random walk {Rn} 
to the chosen direction. It is known that weak con- 
vergence of the distribution of Sn/Bn to the one- 
dimensional standard normal distribution for each cho- 
sen direction guarantees us the weak convergence of the 
distribution of Rn / Bn to the two-dimensional standard 
normal distribution (the theorem of Cramer - Wold [24j). 
Thus, we have the problem of one dimensional conver- 
gence. The distribution of Sj cos tpj is symmetric. If 



b = Esj 



' f(s)ds < 00, 



(10) 



then the distribution of Sn/(ctv N) is approximately 
standard normal due to the classical central limit the- 
orem, where 



cr 2 = Es 2 cos 2 tpx = b Q /2. 



(ID 



Thus, the expectation of \E CZ | after passage of N particles 
is 



En . = rr X I^E = ^L yJnNb (k p a r ), (12) 



where 



b {k p a r )= xK 2 (x)e~ x2/{2k >^ dx. (13) 
Jo 

The function b {k p a r ) is shown in Fig. QJ}. For the dis- 
cussed parameters, b s» 0.4, and the expectation of the 
on-axis longitudinal electric field is 130V/cm. 

Figure [5] illustrates how the actual field differs from the 
expectation value. These are LCODE [25( simulations of 
the plasma response to the full (not half-cut) SPS beam 
[l9| . Five thin curves correspond to five different ini- 
tializations of the random number generator. The thick 
curve is the expectation (|T2|) . The dotted curve is the 
average over 60 different initializations of the random 
number generator. 

The expectation for the off-axis longitudinal electric 
field can be found in a similar way, but with a more 
complicated probability density distribution used instead 
of (Ql . Alternatively, the expectation can be obtained by 
Monte Carlo simulation of the random walk at a reduced 
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FIG. 2. Simulated field amplitudes for various beam ensem- 
bles (thin colored lines), the expectation value (black thick 
smooth line), and the field amplitude averaged over 60 beam 
ensembles (red dotted line) for LCODE simulations of SPS 
beam. 




FIG. 3. The expected value of |_E z (r)| as compared to the 
expectation of |.E S (0)| for oy = kp 1 (thick line) and the cor- 
responding radial profile of the beam density (thin line). 



number of steps (Fig. [3]). Note that the field area is wider 
than the beam itself. 

Since the electric field excited in the plasma is a po- 
tential one, we can find the transverse field component 
from the longitudinal one through the potential $: 



<9$ 

E z = -— = -ikp$, 
oz 



Ej =- 



<9$ 

dr± 



1 dE z 
ik p df± 



(14) 



Typical portraits of the potential $ at the plane (x,£) 
are shown in Fig. [4] These graphs are obtained from 
the random walk model by summing up contributions 
of separate test macro-particles for three different beam 
ensembles. We see that the transverse scale of the poten- 
tial change is shorter than the longitudinal scale. Conse- 
quently, typical transverse fields must be higher than the 
longitudinal ones given by (|12p. 

Calculation of expectation values for transverse fields 
is tricky because of a diverging integral. The model of 
random walk with isotropically distributed steps is still 
applicable, but the absolute value of the step s now de- 
pends of two parameters. For the field component E x we 
have 
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FIG. 4. Maps of the wakefield potential for three different 
beam ensembles. 



where 9 is the polar angle of particle projection onto the 
(x,y) plane. To find the second probability moment b± 
of s, we need the double integration: 



b x = 



'f(s)ds = 



dr 



d6 



s 2 r 



2e 2 k 2 p 



xK 2 (x)e- x2 ^ 2k >^dx = 



4e 2 k% 



(16) 



Here we limit the integration interval by some minimum 
radius r m ; n to avoid divergence. At smaller radii, for- 
mula ([2]) is no longer valid. Several factors could limit 
the applicability of the linear cold fluid model used in 
derivation of @. 

At scales below r% = n p 1 , the plasma cannot be con- 
sidered as a continuous medium. For n p — 7 x 10 14 cm~ 3 , 
r\ w 10 -5 cm. 

The linear theory (23| gives the following perturbation 
of the plasma electron density by a single relativistic pro- 
ton: 



Sn p = k p 6(r± — r&) sinfc p Zf,, 



(17) 



where 5(r± — 7%) is the two-dimensional delta-function. 
If we distribute the same amount of excess charge over 
the area r\ to fulfil the linearity condition 6n p < n p , 
then we find the scale of linearity violation by a point 
charge: T2 — \Jk p Jn p w 3 x 10 _7 cm <C T\. Therefore 
applicability of the linear theory is not a limitation in 
our case. 

Nonzero temperature T e ~ 5 eV of plasma electrons 
modifies the plasma response on the scale of Debye length 



r d w 743, 



/ T e (eV) 
n p (cm -3 ) 



- 6 x 10 5 cm 



(18) 



The modification does not take effect immediately behind 
the beam proton, but after thermally moving electrons 
have time to shift by the distance ~ r^. This occurs at 
the time scale w" 1 , so the Debye length is a limiting scale 
for long proton beams. 

The nonzero plasma temperature results in the nonzero 
group velocity of plasma waves which depends on the 
wavenumber k of the perturbation: 



s(r,6) = 2ek 2 K 1 {k p r)\cos9\, 



(15) 



v g = "&Lu p r 2 d k. 



(19) 
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FIG. 5. a). Growth of the accelerationg field with propagation distance in plasmas. The broken line corresponds to an unseeded 
SPS proton bunch; the solid line shows the self-modulation growth when a 10 MeV electron beam has been used as a precursor, 
b) 2D on-axis cut of the wake field in the unseeded simulation behind the driver at the distance z = 9 m. The field is asymmetric 
due to the competition between hosing and self-modulation; c).2D on-axis cut of the wake field in the seeded simulation behind 
the driver at the distance z = 4 m. The wake is symmetric. 



The sharply localized field spike ([2]) behind the proton 
is produced by short-wavelength wave harmonics, with 
transverse wavenumbers k ~ r^ 1 being responsible for 
the field spike in the area of the transverse size r$. The 
higher k the faster the energy drifts out of the field spike. 
Assuming the time of field evolution er 2 /c, we find the 
minimum scale r% for which the wave has no time to 
drift out of the field spike: 



v g (7 z /c = r 3l r 3 = y 'ik p r 2 d a z « 3 x 10 3 cm « 0.1 k p x . 

(20) 

This is the scale we take for r m ; n . The function bx(k p a r ) 
for k p r m - m = 0.1 is shown in Fig. QJ>. For a r = k~ x , 

b\ ~ 0.85, and the expectation of the on-axis transverse 
electric field component is 



E„ 



irNbi(k p a r ) « 200V/cm. (21) 



An uncertainty in determination of r m ; n has a little effect 
on the field expectation because of the weak (logarithmic) 
dependence of b\ on r min . Factor of two smaller r m ; n 
results in 0.5 In 2 (or 40%) increase of b\ (Fig. [TJj) and 
20% increase of E ax . 

We have simulated the wake field by the SPS beam 
using the hybrid 3D particle- in-cell code VLPL [H,^. 
We simulated both the "unseeded" case, when the wake 
field grows from the shot noise of the SPS bunch, and the 
case when the proton bunch self-modulation is seeded by 
a precursor electron bunch. The major difficulty is to 
simulate the shot noise properly. The number of numeri- 
cal macroparticles is several orders of magnitude smaller 
than the 1.15 ■ 10 11 protons in the SPS bunch. Thus, 
the noise generated by randomly seeded macroparticles 
would be significantly higher than the natural expecta- 




tion (|12p . For this reason, we initialize all the numer- 
ical beam macroparticles at the centra of the grid cells 
and distribute their transverse momenta regularly and 
symmetric according to the bunch divergence. The lon- 
gitudinal momenta have a random Gaussian distribution 
centered at 450 GeV/c and the spread corresponding to 
the given longitudinal omittance. A bunch initialized in 
this way generates an extremely low numerical noise field, 
a couple orders of magnitude below the expectation (fl2|) . 

To adjust the noise to the natural level (fP2")) . we intro- 
duce random displacement e to the regular positions of 
macroparticles in the configuration space: 



(22) 



Here r is the space coordinate, h is the cell size, Np is the 
number of numerical macro-particles which substitute for 
Np real beam protons in the cell volume, and np is the 
local beam density. The random displacements (f22|) lead 
to the same level of numerical beam density fluctuations 
as the y/Np fluctuations of the real beam proton number 
within the cell volume. 

Fig. [5] shows results of the 3D simulation. The frame 
Fig. [SJa) gives the wake field growth with propagation 
distance. The broken line corresponds to the simulation 
without any external seeding. The proton bunch self- 
modulation started to grow from the natural beam den- 
sity fluctuations due to the shot noise. We see that it 
takes 10 m of plasma for the wake to reach its maximum 
at about 80 MV/m. At the distance of 10 m, the proton 
bunch has already significantly diffracted due to its trans- 
verse emittance. The generated wake field behind the 
proton bunch is shown in Fig. [SJb) . We zoom here at the 
position 2er z behind the middle of the SPS bunch. One 
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sees that the wake is not very symmetric and is tilted. 
This tilt is due to the presence of the hosing instability 
that competes with the axisymmetric self-modulational 
mode The solid red line in Fig. 03a) shows the wake 
field growth when the self-modulation is seeded by a pre- 
cursor electron bunch. We used a 10 MeV electron bunch 
of 1 ps duration and the current of 100 A as the pre- 
cursor. It is seen that the field fluctuates strongly at 
the first meter of propagation in plasma. This is due 
to strong self-focusing/defocusing of the electron bunch. 
After the first meter, the self-modulation of the proton 
bunch sets in and the field starts to grow exponentially. 
The maximum field of some 0.6 GV/m is reached after 
4 meters propagation distance. The field snapshot at 



this distance is given in Fig. Elc). Apparently, the field 
is symmetric so that the self-modulation mode wins over 
the hosing. Our simulations also show that the growth of 
self-modulation is not very sensitive to the electron bunch 
parameters such as the energy and the current. It is im- 
portant that the initial wake generated by the precursor 
is significantly larger than the noise field of the proton 
bunch and that the saturation of the self-modulation is 
reached before the proton bunch diffracts away due to its 
transverse emittance. 
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